Muestreo de Fósforo y Nitrato

Campaña de muestreo

Las muestras de Fósforo y Nitrato se hace mediante campañas esporádicas con frecuencia de ~2 meses. En dichas campañas no se realizan aforos. A continuación se presentan los valores muestreados en el período 2004-2016 así como el caudal ajustado que de cada día (ver: tratamiento de caudales).

library(kableExtra)
qlt_41 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/data/rch41_calidad.RDS")
qlt_41$qobs = qlt_41$qobs*1.12
kable(qlt_41[,1:4], booktabs = TRUE,
             caption = "Muestras de Fósforo y Nitrato en Paso Pache",
      col.names = c("Fecha", "Nitrato (mg/l)","Fósforo (mg/l)", "Caudal (m3/s)")) %>%
  kable_styling(c("striped", "bordered")) %>%
  scroll_box(width = "600px", height = "200px")

Muestras de Fósforo y Nitrato en Paso Pache
	Fecha	Nitrato (mg/l)	Fósforo (mg/l)	Caudal (m3/s)
1	2004-12-21	0.050	0.1600	5.83856
2	2005-01-25	0.230	0.1500	2.90640
3	2005-02-22	0.340	0.2100	6.28880
4	2005-04-19	0.150	0.1800	280.94304
5	2005-08-02	0.630	0.1200	37.36880
6	2005-11-15	0.050	0.4900	13.58112
7	2006-03-21	0.530	0.1500	38.51456
8	2006-07-04	0.320	0.0700	49.92736
9	2006-09-26	0.130	0.1100	23.18064
10	2007-01-23	0.510	0.0900	0.89824
11	2007-06-05	0.550	0.4400	15.77968
12	2007-09-11	0.570	0.1800	40.75680
13	2007-11-20	0.050	0.1100	10.61648
14	2008-01-22	0.850	0.2400	2.28032
15	2008-05-20	0.050	0.1200	3.38352
17	2009-03-17	0.790	0.2200	79.80336
18	2009-04-28	0.060	0.0800	1.36528
19	2009-08-26	0.400	0.1900	33.30208
20	2010-05-19	0.070	0.1600	4.41392
21	2010-07-20	0.220	0.3600	490.98000
22	2010-09-28	0.320	0.1400	24.90320
23	2010-11-16	0.470	0.1400	11.40832
24	2011-01-20	0.014	0.7030	1.97568
25	2011-03-17	0.090	0.2580	2.17056
26	2011-05-05	0.035	0.1390	4.50688
27	2011-07-21	0.330	0.2430	98.70336
28	2012-01-04	0.014	0.3090	2.21536
29	2012-03-08	0.130	0.2250	2.84144
30	2012-05-10	0.015	0.1580	2.71936
34	2013-01-10	0.230	0.2400	20.43440
35	2013-03-14	0.260	0.3200	19.37488
36	2013-07-25	0.430	0.2790	20.93280
37	2014-01-23	0.350	0.2220	17.04752
38	2014-03-20	0.290	0.2000	18.72416
39	2014-05-08	1.000	0.2070	22.21520
40	2014-07-10	0.360	0.1880	116.85856
41	2014-09-25	0.370	0.1920	24.79456
42	2014-11-13	0.020	0.0973	24.13824
43	2015-01-22	0.380	0.2000	4.70848
44	2015-03-19	0.016	0.2920	2.16160
46	2015-07-23	0.020	0.1290	2.17840
47	2015-09-24	0.200	0.1700	6.73680
48	2015-11-19	0.055	0.1450	4.25040

qlt_23 <- readRDS("~/R_projects/SWAT_calibration_ETmodis/data/rch23_calidad.RDS")
qlt_23$qobs = qlt_23$qobs*1.65
kable(qlt_23[,1:4], booktabs = TRUE,
             caption = "Muestras de Fósforo y Nitrato en Paso de los Troncos",
      col.names = c("Fecha", "Nitrato (mg/l)","Fósforo (mg/l)", "Caudal (m3/s)")) %>%
  kable_styling(c("striped", "bordered")) %>%
  scroll_box(width = "600px", height = "200px")

Muestras de Fósforo y Nitrato en Paso de los Troncos
	Fecha	Nitrato (mg/l)	Fósforo (mg/l)	Caudal (m3/s)
1	2011-01-20	0.014	0.5260	0.11880
5	2012-01-04	0.014	0.1300	0.11880
6	2012-03-08	0.300	0.2090	6.53235
7	2012-05-10	0.080	0.0752	0.68310
8	2012-06-28	0.100	0.0905	2.29515
9	2012-09-20	0.340	0.0975	47.77905
10	2012-11-08	0.120	0.1140	8.29620
11	2013-01-10	0.190	0.0838	1.43550
12	2013-03-14	0.140	0.0745	1.85790
13	2013-07-25	0.310	0.0869	5.28000
14	2014-01-23	0.150	0.0892	6.35910
15	2014-03-20	0.110	0.1000	3.15975
16	2014-05-08	0.160	0.0748	5.35260
17	2014-07-10	0.450	0.0566	18.97335
19	2014-09-25	0.130	0.0906	7.84245
21	2014-11-13	0.180	0.1300	8.65260
22	2015-01-22	0.190	0.0780	1.82985
23	2015-03-19	0.020	0.0830	0.11880
25	2015-07-23	0.020	0.0450	0.41085
26	2015-09-24	0.110	0.0830	4.80975
27	2015-11-19	0.094	0.0500	0.87450
28	2016-01-21	0.130	0.1270	0.11880
29	2016-04-07	0.110	0.1010	2.46675
30	2016-06-30	0.090	0.0700	32.92080
31	2016-10-13	0.083	0.1570	2.32320

Relación de Fósforo - Caudal

La siguiente figura presenta la relación entre el caudal unitario (lps/km²) y la producción unitaria de fósforo (mg/s/km²). Ha sido obtenida a partir de los muestreos de calidad de agua en Paso Pache y Paso de los Troncos. Los valores de caudal has sido tomados de las correspondientes estaciones de caudal.

library(plotly)
area_41 = 4896
area_36 = 2744
area_23 = 687
qlt_41$qobs = 1000*qlt_41$qobs/area_41
qlt_23$qobs = 1000*qlt_23$qobs/area_23
qlt_41$NO3_mgxs = qlt_41$NO3_mgxs/area_41
qlt_23$NO3_mgxs = qlt_23$NO3_mgxs/area_23
qlt_41$PT_mgxs = qlt_41$PT_mgxs/area_41
qlt_23$PT_mgxs = qlt_23$PT_mgxs/area_23
fit41_PT = lm(log(qlt_41$PT_mgxs)~log(qlt_41$qobs))
fit41_NO3 = lm(log(qlt_41$NO3_mgxs)~log(qlt_41$qobs))
fit23_PT = lm(log(qlt_23$PT_mgxs)~log(qlt_23$qobs))
fit23_NO3 = lm(log(qlt_23$NO3_mgxs)~log(qlt_23$qobs))
fit41_PT = round(c(exp(fit41_PT$coefficients[1]), fit41_PT$coefficients[2]),2)
fit41_NO3 = round(c(exp(fit41_NO3$coefficients[1]), fit41_NO3$coefficients[2]),2)
fit23_PT = round(c(exp(fit23_PT$coefficients[1]), fit23_PT$coefficients[2]),2)
fit23_NO3 = round(c(exp(fit23_NO3$coefficients[1]), fit23_NO3$coefficients[2]),2)
qlt_41$fit_PT = fit41_PT[1]*qlt_41$qobs^fit41_PT[2]
qlt_41$fit_NO3 = fit41_NO3[1]*qlt_41$qobs^fit41_NO3[2]
qlt_23$fit_PT = fit23_PT[1]*qlt_23$qobs^fit23_PT[2]
qlt_23$fit_NO3 = fit23_NO3[1]*qlt_23$qobs^fit23_NO3[2]

PT <- plot_ly(data = qlt_41, x = ~qobs, y = ~PT_mgxs, type="scatter",
              mode   = 'markers',colors="blue", width=2, name="Paso Pache")
PT <- PT %>% add_lines(data=qlt_41,x = ~qobs, y = ~fit_PT,
                         type = 'scatter',
                         mode = 'line', 
                         line = list(color = 'blue', width = 2),
                       name=paste0("PT = ",fit41_PT[1],"Q^",fit41_PT[2]))
PT <- PT %>% add_markers(data=qlt_23,x = ~qobs, y = ~PT_mgxs,
                       type = 'scatter',
                       mode = 'markers', 
                       marker = list(color = 'red', width = 2),
                       name="Paso de los Troncos")
PT <- PT %>% add_lines(data=qlt_23,x = ~qobs, y = ~fit_PT,
                       type = 'scatter',
                       mode = 'line', 
                       line = list(color = 'red', width = 2),
                       name=paste0("PT = ",fit23_PT[1],"Q^",fit23_PT[2]))

PT <- layout(PT, xaxis = list(title="Caudal (lps/km2)",range=c(-1,2),type = "log"),
             yaxis = list(title="PT (mg/s/km2)", range=c(-2,1),type = "log"),
             legend = list(x = 0.01, y = 1),
             width=500,
             height=500

)

PT

Relación Nitrato - Caudal

La siguiente figura presenta la relación entre el caudal unitario (lps/km²) y la producción unitaria de Nitrato (mg/s/km²). Ha sido obtenida a partir de los muestreos de calidad de agua en Paso Pache y Paso de los Troncos. Los valores de caudal has sido tomados de las correspondientes estaciones de caudal.

NO3 <- plot_ly(data = qlt_41, x = ~qobs, y = ~NO3_mgxs, type="scatter",
              mode   = 'markers',colors="blue", width=2, name="Paso Pache")
NO3 <- NO3 %>% add_lines(data=qlt_41,x = ~qobs, y = ~fit_NO3,
                         type = 'scatter',
                         mode = 'line', 
                         line = list(color = 'blue', width = 2),
                       name=paste0("NO3 = ",fit41_NO3[1],"Q^",fit41_NO3[2]))
NO3 <- NO3 %>% add_markers(data=qlt_23,x = ~qobs, y = ~NO3_mgxs,
                       type = 'scatter',
                       mode = 'markers', 
                       marker = list(color = 'red', width = 2),
                       name="Paso de los Troncos")
NO3 <- NO3 %>% add_lines(data=qlt_23,x = ~qobs, y = ~fit_NO3,
                       type = 'scatter',
                       mode = 'line', 
                       line = list(color = 'red', width = 2),
                       name=paste0("NO3 = ",fit23_NO3[1],"Q^",fit23_NO3[2]))

NO3 <- layout(NO3, xaxis = list(title="Caudal (lps/km2)",range=c(-1,2),type = "log"),
               yaxis = list(title="NO3 (mg/s/km2)", range=c(-2,1),type = "log"),
             legend = list(x = 0.01, y = 1),
             width=500,
             height=500
)

NO3